import pymap3d as pm

# The local coordinate origin (Zermatt, Switzerland)
lat0 = 60 # deg
lon0 = 120  # deg
h0 = 0     # meters

# The point of interest
lat = 31.0  # deg
lon = 121.0   # deg
h = 10      # meters

print(pm.geodetic2enu(lat, lon, h, lat0, lon0, h0))